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We investigate magento-oscillations in the specific heat of non-relativistic fermions with a Fermi 
surface minimally coupled to a fluctuating U(l) gauge field. Our study is motivated by the recent 
observation of quantum oscillations in the underdoped cuprates, and by theoretical models of pocket 
Fermi surfaces realizing a non-Fermi liquid 'algebraic charge liquid'. Our main result is the compu- 
tation of the order 1/N correction to the Lifshitz-Kosevich expression for the oscillation amplitude 
in the dirty limit in a model with N species of fermions. 



I. INTRODUCTION 

The problem of two dimensional non-relativistic 
fermions coupled to a gauge field has been intensively 
studied over the last couple of years in a variety of dif- 
ferent contexts. It appears as a low energy description 
of different models of strongly correlated electronic sys- 
tems, such as electrons in the fractional quantum Hall 
regime and in theories of non-Fermi liquid phases for the 
underdoped cuprates. 

Our study is motivated by the remarkable recent 
observation a 1 ' 2 ' 3 ' i 5 ' 6 of quantum oscillations in the un- 
derdoped cuprate superconductors at high magnetic 
fields. So far, these observations have been consistently 
interpreted using Fermi liquid models of the formation 
of Fermi pockets. The Fermi liquid theory of quan- 
tum oscillations^ yields periodic behavior as a function of 
1/-B, where B is the applied magnetic field, with an am- 
plitude given by the Lifshitz-Kosevich (LK) prefactor 5 . 

However, as the precision and range of the observa- 
tions increase, it would be useful to have theoretical pre- 
dictions for other candidate metallic ground states of 
the underdoped cuprates. To this end, we will examine 
the amplitude of the quantum oscillations in 'algebraic 
charge liquids' (ACL) ^ 10 ' 11 . The charged excitations in 
these states are described by Fermi surfaces of spinless 
electrons coupled to an emergent U(l) gauge field, as we 
will review in more detail in Section [Vj We will find that 
these systems also exhibit oscillations which are periodic 
in 1/-B, with a prefactor with small but detectable devi- 
ations from the LK theory. 

Ref. 12 has also addressed this problem in the clean 
limit. However, they do not include the oscillatory terms 
in the gauge field propagator, which are responsible for 
the main effects we describe below. 

It is also possible that our results have implications 
for quantum Hall systems, where similar theories apply 
to compressible states at even denominator fillings. How- 
ever, we will not explore this connection further here. 

Our main analysis is a computation of the free en- 



ergy of a system of N species of fermions coupled to a 
fluctuating, emergent U(l) gauge field. We compute the 
de Haas-van Alphen oscillations in the presence of an 
applied magnetic field in the dirty limit to order 1/N. 
The N = oo theory co-incides with the LK result. Our 
main results for the order 1/N corrections are shown in 
Eq. dUD, Eq. ([H, and Eq. flU}. Most interestingly, we 
find a qualitative difference in the behavior of the os- 
cillations as a function of T/u> c , i.e. temperature over 
the cyclotron frequency, compared to Fermi liquid the- 
or y 13 ' 14 i 15 . This comparison is plotted in Fig. [TJ Addi- 
tionally, for the specific heat, cy, we find that the gauge 
field correction to the oscillatory term has a temperature 
(T) dependence ~ Tln(l/T) which differs from the ~ T 
dependence in the LK term, and so may be experimen- 
tally detectable in recent and future experiments 16 . 

The organization of the paper is as follows. In Sec- 
tion HD we introduce the generic model of fermionic de- 
grees of freedom in the presence of weak scalar disorder 
coupled to an internal gauge degree of freedom. In Sec- 
tion UTS] we discuss the Landau level structure, and the 
role of disorder together with its treatment on a phe- 
nomenological level. We conclude this section by deriv- 
ing the effective theory for the gauge field in the presence 
of disorder and Landau levels in Section Hi Bl 



In Section IIIII we present a formalism to extract the 
oscillatory part of the gauge field contribution to order 
1/N. To this end we have to calculate oscillatory ther- 
modynamic and transport quantities, part of which is 
outlined in Appendix [XJ In Section IIVI we show how to 
calculate the specific heat of the entire coupled system to 
order 1/N. In Section ITVAl we review the derivation of 
the specific heat and the associated de Haas-van Alphen 
oscillations of a disordered gas of electrons subject to a 
magnetic field. In Section TlVBI we calculate the specific 
heat of the gauge field. In a first step (Section IIVB ip 
we derive the specific heat of the non-oscillatory part, 
whereas in Section llV B 2l we calculate the oscillatory cor- 
rection. In Section [V] we explain the meaning of our re- 
sults in the context of a very recent spin density wave 
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(SDW)+ gauge field description of the cuprates. Finally, 
we summarize in Section [VII 



II. MODEL 



magneto-oscillations, which is what we are interested in, 
is characterized by uj c t <C 1. In this regime the density 
of states, for instance, has a constant part with a smooth 
oscillatory part on top of it, see Eq. (fT5|) . 



The effective model which forms the basis for this pa- 
per was introduced in different contexts and its main 
characteristics are reviewed here. We consider a con- 
ventional Fermi gas with quadratic dispersion, which is 
minimally coupled to an internal U{1) gauge field. The 
Lagrangian for this system generically reads 



(V - iKf 
2m 



where A denotes the internal U(X) gauge field. 



(1) 



Disorder broadened Landau levels in the dirty 
limit 



In our problem we consider the above system subject 
to an external perpendicular magnetic field B = V x a, 
which implies the fermions also couple minimally to the 
external gauge field a and causes the electronic system to 
organize its spectrum into Landau levels. In the absence 
of the internal U(l) gauge field A we can diagonalize the 
electronic part and cast it as 



f d T 



f 



(2) 



where uj c = eB/m is the cyclotron frequency, n the Lan- 
dau level, and the degeneracy of the Landau levels is 
given by 



1 



2n 



UJ C = V UJ C 



(3) 



where Vq = is the density of states of the two dimen- 
sional electrons at the Fermi level. 

In addition, we want to consider dilute disorder which 
couples to the electromagnetic charge but not to the spin 
sector. A consistent treatment of such a system was 
presented in the context of the two dimensional disor- 
dered electron gas in a perpendicular magnetic field by 
Andoi£, who incorporated this in the framework of the 
self-consistent Born approximation. On a phenomeno- 
logical level in the oscillatory regime this can be done 
by introducing a finite lifetime for the electronic degrees 
of freedom. The one particle retarded Green's function 
reads 



G f {uj,n) 



1 



■j c (n + \) + 



(4) 



The major effect of the disorder is thus to broaden the 
Landau levels into a Lorentzian shape. The regime of 



B. Effective action for the gauge field in the 
diffusive limit 

In order to calculate the specific heat in a 1/N- 
expansion, we need to derive the photon propagator of 
the internal U{1) gauge field in the presence of disorder 
and Landau levels. Schematically, our derivation goes 
along the following line: the fcrmionic action reads 

S = drtfxf [id T +fi- e(-iV - ea — A) + iA T ] f . 
Jo 

(5) 

We assume a parabolic dispersion of the fermions, i. e. 
e (k) = In a next step we integrate out the fermions, 
leading to 



<S c ff = — tr In [id r + fi — e(— iV — ea — A) 



(6) 



Expanding the above expression to second order in A 
using — iV — ea = II we obtain the effective gauge field 
propagator from the polarization operator. In principle, 
one could derive the effective polarization operator in the 
presence of disorder in the basis of the Landau levels^. 
However, here we choose a different route. 

In general, the polarization operator is obtained by 
expanding to second order in the internal gauge field and 
performing the functional derivative according to 



->cir 



<5A l (q, ii/)5A J (— q, —iv) 



= lP(q> n ) 



(7) 



Following Halperin et alJ& we keep the following two- 
component form of the photon-propagator 



D-^q.w) =nq,w) = 



IT 00 



q_Yl x v 
f\yy 



(8) 



Subsequently, we concentrate on the low energy form of 
the photon propagator in the diffusive regime. The form 
of the photon propagator is highly constrained by con- 
servation laws and generically reads 



n 00 

fpv 

Il xy 



Dq 2 



Dq 2 — iu> 
iuia vy + q 2 x , 

i,.,rr x V 



(9) 



The diffusive limit of the system is accounted for by re- 



placing II 00 = vq 



Dq 



D > q - which on a formal level is 
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achieved including impurity ladders in the vertex func- 
tion. The dependence upon the magnetic field enters 
through (j xx , cr xy , v, and \- D denotes the diffusion con- 

2 

stant and is set by D — ^pr in two dimensions. The free 
energy due to the gauge field is readily calculated (see 
Halperin et al*2-) 



f A = 



A d 2 q f duj , ^ Im dct £> _:L (q,cj) 



(2 



~To / — nb(cj)arctan 

TT) Z J 2lT 



Re det D^^uj) 
(10) 



The upper cutoff A is set by roughly twice the Fermi 
momentum, thus A w 2fci?, which is the upper bound 
for the existence of low-energy excitations. Eq. (jTUJ) is 
the central expression which will allow to calculate the 
contribution of the internal gauge held with and without 
applied external magnetic held. 



III. EXTRACTING THE OSCILLATORY PART 
OF THE SPECIFIC HEAT OF THE GAUGE 
FIELD 

In the following we decompose the inverse photon prop- 
agator (Eq. ([8])) into two parts, one containing the non- 
oscillatory contributions, called D and another part con- 
taining the oscillatory contributions, called D osc . The 
quantities a xx , a xy , x, D, and v naturally decompose 
into a non-oscillatory part and an oscillatory part 



@xx 

cr xy 



i/ ■ 

_0 



xx 




x = x° + x osc 



(11) 

where all oscillatory contributions are expressed as a 
power series in exp ^-^77^, which on a formal level is 
obtained through a Poisson summation duality. As long 
as oj c t 1 it suffices to retain the first moment in this 
power series to isolate the leading oscillatory contribu- 
tion. The results for the oscillatory components are given 
in the subsequent section (Sec. IIII A]) . 

Using the property from standard perturbation theory 

det I)- 1 = dct D Q Met (1 + D ■ £>~c) 

« detD^il + trDo-D-l) (12) 

we can formulate the contribution of the gauge held to 
the free energy as 



/ = 



O 



A d 2 q 


f dto 


(2tt)2 J 


' 2T 


A d 2 q 


f dui 


(27T) 2 J 


' 27 


A d 2 q 


f dto 


(2tt)2 J 


' 27 


(.-*) 





rib(u)) arctan 



nb(u}) arctan 



Im dct D _1 (q, 
Re det D- X (q,a;) 
Im det l)(^ 1 (q, uj) 
Re dct D^^uj) 



(13) 



It is important to note that the second term is now pro- 
portional to exp {^-^j-^j and thus small compared to the 

first term. The first term in the above expression has 
been analyzed by Halperin et alJ^, and is known to yield 
a contribution to the free energy of the type T In T, which 
is reviewed later. The additional factor of ln(I/T) is the 
manifestation of the well-known Altshuler-Aronov^ cor- 
rection. 

In the following we retain the leading order temper- 
ature dependence of the density-density and current- 
current response in fi/T, however we will allow for ar- 
bitrary T/uj c . 



A. Oscillatory thermodynamic and transport input 
quantities 

The DOS can be calculated in a way analogous to the 
grand potential (see Appendix |A"]). At the Fermi level, in 
the regime r/i 3> 1, the expression for the DOS reads 



1=1 



,( M ) = V0 [l + 2 y £(-l) l co B ^e-&) , (14) 



whose leading oscillatory behavior in the limit lo c t <C I 
reads 



27T/1 

s» vq I 1 — 2 cos e 



(15) 



We can calculate the longitudinal conductivity accord- 
ingly 17 , yielding 



OO 



a» = .o^^[ 1 + 2^-1)' cos 27 ^ 6 '^ 



00 



' 1 + (UJ C T) 2 
1 

1 + (uJcT) 2 



1 = 1 



uj r sinh A 



' 2-Kfj, — 2_ Ai 
1 — 2 cos e w r 



sinh Ai J 



(16) 



with go = — and A; = 

U TCI 1 UJ C 

Using the well known Einstein relation for diffusive sys- 
tems 



a = vD 



(17) 



we can determine the oscillatory part of the diffusion con- 
stant D. It turns out that to leading order we have 

D = D ^l + 2cos^V^ (1-C(7> C ))) , (18) 

where we introduced the function 

2it 2 T/uj c 



C(7>c) 



(19) 



sinh 2tt 2 T/uj c 

for notational convenience. We observe that for T <C w c 
D = Do to all orders in exp ■ Furthermore, it is 

straightforward to show that 

&xy = —U c T(Jxx , (20) 
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to all orders in exp ^— jrf J ■ 

The diamagnetic susceptibility can be obtained from 
the grand potential (Eq. (|A14|) 1 by 



d 2 n 
ob 2 



(21) 



Using Eq. (|A14|) in the limit lo c t < 1, < 1, and 
/jr > 1 we obtain 

1 



X. 



24-7TTO 



[I 2 2TTfi 



1 + 24^r cos 

[l 2 2TT/J, 



Xo 1 + 24^ cos 



C(T,w c )) . (22) 



where 
a(T) 



<9 2 

4tt 6 T 
Wc sinh 3 (^) 

4^4 sinh C i^fT 



sinh 



3 + cosh 



/ 47T 2 T 



(29) 



This function is dimcnsionless and compared to the os- 
cillatory contribution from the gauge field in Fig. [TJ 
In the limit — <C 1 this reduces to the well-known 



IV. FREE ENERGY AND SPECIFIC HEAT 



^ 2 ^o . 
c v = ~^—T | 1 



2 cos e 



(30) 



In the spirit of the large- TV treatment we can expand 
the free energy T of the system to order 1/N, which 
yields the following result 



T = Nf f + f J 



(23) 



where J -' is the free energy of the non-interacting 
fermionic system and f A denotes the free energy asso- 
ciated with the fluctuations of the emergent gauge field. 

The specific heat of this expression can be obtained by 
the well-known formula cy — — T^§. This implies that 
the specific heat decomposes into two parts 

cy = Nc f v + 4 , (24) 

which are analyzed independently. 

A. Specific heat of the electrons 

The free energy is related to the grand potential via 
Legendre transform according to 

ff = f iN + n, (25) 

which allows to obtain the specific heat via 

03/0/ 



cv 



-T- 



dT 2 



(26) 



The dependence of the chemical potential upon the mag- 
netic field is subdominant, thus 



cv = —T 



,9Hl 
df 2 



(27) 



with the oscillatory contribution obtained from Eq. (|A14D 



, . v 2tth _ 
oiyl ) — ^2 cos e 

a{T)— cos— f=-e" 

47T" 3 Ul c 



(28) 



per spin species. 



B. Specific heat of the gauge field 

In the following we calculate the non-oscillatory and 
oscillatory contributions of the gauge field to the specific 
heat separately. 



1. Non-oscillatory contribution 

We start with a calculation of the non-oscillatory part 
of the specific heat of the gauge field. In order to do so 
we analyze 



d 2 q 



dw det D^H^uSf 

— rihiw) arctan . (31 

2n detD^( q ,cj)' 



We can rewrite the contribution to the specific heat as 
T r A/VT 



dk 



16tt 2 

k(2x-x 2 coth(|)) 
sinh(a:/2) 2 



dx 



arctan 



cxk 



ak 4 + bx 2 



(32) 



with 



a 
b 
c 



D XQ (l + 48Tr 2 (u c T) 2 (a v y ) 2 ) 



'0 ) 



Xo 



Do? 



(33) 



for lo c t <C 1 and /it 3> 1. In order to analyze the asymp- 
totic behavior of the above ir-integral, we start noting 
that the x-integration is cut off by the factor 



Sinh 2 ■ 
1 



the order of x = 10. If we consider the factor ak i+ bx i , we 

t) 1 ^ 4 it becomes -4r, whereas 



know that for k 3> v 10 I 
for k <C 4 it becomes -i? 



Only in the former case will 
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the integral contribute a logarithmic dependence upon 
temperature, hence 



-v 



cT 
6a 



a/Vt 

Ml) 



dk 

1/4 k 



cT , AV/ 2 
In ■ 



12a 1Q6V2T 



We furthermore introduce the short form 

K=l+48(fiT) 2 (uj c T) 2 , 

which is a number of order 1. This leaves us with 
as the non-oscillatory contribution. 



(34) 



(35) 



(36) 



2. Oscillatory contribution of the gauge field 

In this section we analyze the oscillatory contribution 
of the gauge field. 

Following the prescription given in Eq. (|12p we decom- 
pose the gauge-field propagator according to 



\ -iq<r xy iujCT yv + X q 2 J 

( Dnq 2 ■ xy 

= v o Do< p_ iu) iq<r u 

I ■ xy ■ yy , 9 

\ -iqo- y iLua^+xoq 



Dq 2 



_|_ I y Dq 2 —iu) 



osc 



(37) 



Using the oscillatory expressions derived in Section IIII Al 
we identify 



27T/i 



D^l ~ — 2 cos e "c Dr. + cos e "c 



27T/i 



<5£f (i - C(i> )) 



24gxog 2 C(r,^c) / ' 

(38) 



where ^ <§; 1 was used. Calculating Im tr £>o ■ D os \ 
we realize that the first term on the right hand side of 
Eq. (|38j) does not contribute an imaginary part, thus 

- , ti 2 27T/1 

Im tr D • -Dose = -24^-tt cos — -e ^ x 



I , , r , C( i 2uj 



Ao 4 + £w 2 D 2 g 4 + u 2 



(39) 



Looking at Eq. (|39|) it is interesting to note that there are 
two contributions. The first contribution survives in the 
limit T/u>c — > 0, whereas the second part goes to zero. 

Additionally, the second term is down by a factor — §■ and 
is thus parametrically small compared to the first term. 
Consequently, we will discard the second contribution for 
our analysis. The remaining constants read 

A = D 2 (m 2 -^oXo) 2 , 

b = ((C) 2 + K ra f) 2 «W) 4 , 

,2 



(40) 



The constants of the disordered electronic system in- 
volved in the above expression assume the well-known 
values 

m „2 



,yy _ 



1 



1 



to 1 + (lu c t) 2 1 + (w c r) 2 

UJ C T 



CT » (TO , 



1 + (w c t) : 



Ma 

,,2 



m 
2^ 



1 



2tt 



Xo 

2 



247TTO ' 
2 _ 27T(To 



'0 > 



77 J 



(41) 



We now calculate the temperature dependent part of the 
free energy (see Appendix [BJ . This part reads 



H 2 C 



/osc — 2 A C0S e 



H 2 C 27T/J 

cos e -orC(T,a; c )T 



(42) 



with £ w 0.352822. We introduce the dimensionless func- 
tions 

= -^T 2 C(T,^ C ), 

7 (T) = -C(T, Wc ) + ^T 2 C(T )Wc ) + ^(T) ) (43) 
which allows to express the specific heat as 



/lose 
C V 



27r/i 24 (/xr) J 

cos e »c 2^ — toJ m 

w c (u} c t) k 2 

2ttu zr 24 (ht) 3 , , 



'A A 2 
B T 



(3(T) 



(44) 



with k as defined in Eq. ([35]). If we compare the oscil- 
latory contribution of the LK type (Eq. ([28]) ) with the 
1/AT corrections we realize that their ratio behaves as 



„Aosc 



(T) 96nH^r^ T )^{^^)+l(T) 



c^ c (T) N{oj c t) 2 k 2 



a(T) 



(45) 



6 



0.6 0.8 1.0 



FIG. 1: This plot shows the function a(T) and 50/3(T), where 
T is measured in units of uj c and |£ = 100 in units of u> c . We 
magnified the function /3(T) by the arbitrary factor 50 for a 
better comparison of the two curves. We see that both curves 
saturate for very low temperatures. 



For low temperatures the logarithmically diverging part 
dominates the 1/N correction, which is why we plot the 
functions a(T) and f3(T) only. For a reasonable qualita- 
tive comparison of the functional forms of the two con- 
tributions we take the dimensionless prefactor to be of 
the order 50. Assuming this number, Fig. [T] shows a 
comparison of the temperature dependent dimensionless 
functions a(T) and 0(T). However, one still has to keep 
in mind that the 1/iV-correction has an additional log- 
arithmic temperature dependence, which was scaled out 
for the comparison. 



V. APPLICATION TO THE UNDERDOPED 
CUPRATE SUPERCONDUCTORS 

Recent experiments performed in the underdoped 
regime of the cuprate superconductors show great 
promise to shed light onto these still mysterious ma- 
terials. In our discussion we focus on quantum 
oscillations measurements in the underdoped region 
of YBa 2 Cu 3 07-5 (YBCO)!^^. i n this context, 
LeBoeuf et observed signatures reminiscent of pock- 
ets of carriers of charge — e (in contrast to hole- like 
charge-carriers) . 

Here, we will discuss these experiments using a spe- 
cific theoretical model ^ 10 ' 11 ! 21 ' 22 for the interplay be- 
tween spin density wave (SDW) and superconductivity 
in the underdoped regime. We investigate the magneto 
oscillations for magnetic fields greater than H C 2, i.e. in 
the absence of superconductivity. In contrast to ear- 
lier works 1 ^ which investigated Shubnikov-de-Haas oscil- 
lations, we emphasis the role of de-Haas-van- Alphen os- 
cillations, i.e. oscillation in thermodynamic quantities. 

The theory has two non-superconducting metallic 
phases. One has long-range SDW order, and so is a con- 
ventional Fermi liquid at low enough temperatures: here 
the magneto oscillations will be given by the LK theory. 



The other metallic phase has no SDW order, but retains 
aspects of the Fermi pocket structure of the SDW-ordered 
phase: this is the 'algebraic charge liquid'. The ACL has 
an emergent gapless U(l) photon which will lead to cor- 
rections to the LK theory, as we have discussed above. 
The photon acquires a Higgs mass, A.af, across the tran- 
sition from the ACL to the SDW phase, and so its effects 
are quenched in the SDW phase. 

The specific model has the Lagrangian 



C = C z +C a + C 



f 



(46) 



The first term describes the magnetic degrees of free- 
dom. Conventionally, the slow magnetic degrees of free- 
dom arc expressed in the framework of 0(3) non-linear 
sigma model. Here, however, we map the spin density 
wave (SDW) order parameter to bosonic degrees of 
freedom z, which carry spin S = i, via 



2^ z a?apZp 
a, f3 



(47) 



and the effective Lagrangian assumes the form of the so- 
called CP 1 -model 



N 



(\(dr - iA T )z a \ 2 + v 2 \(X7- iA)z a \ 



,a=l 
N 



+ ip J2(K\ 2 -n) 



a=l 



(48) 



where A T and A denote an internal U(l) gauge field 
emerging from the redundant parametrization of the 
SDW order parameter shown in Eq. (|47|) . Interestingly, 
the bosonic spinons z locally determine the spin axis of 
the physical electrons, which leads to a fractionalization 
of the spin and the charge degree of freedom. Conse- 
quently, the effective charge carriers also couple to the 
internal gauge field. It was known for a long time that 
the existence of spin density wave order, i.e. (z) ^ is 
responsible for a Fermi surface reconstruction 2 ^ 24 . In a 
very simple approximation we can take this fact into ac- 
count by introducing g + and g„ for the electronic pockets 
sitting at the antinodal points ((ir, 0))and /+2, f—i, 
and f-2 for the hole degrees of freedom sitting at the 
nodal points ((tt/2,it/2)): 

fn „ (V-iA) 2 \ 
£ g = 9+ [Or - iA T — fl\g + 

/_ .. (V + zA) 2 \ 
+ 9- [d T + iA T - — - fij g- 



(49) 



and 



= ( «r - nAr - ^ ^ ,> j 



2m 2 



qa i 



(50) 
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where q = ± and a = 1,2. As explained before, all the 
fermionic degrees are coupled to the internal £7(1) gauge 
field, but carry different charges under the transforma- 
tion. 

The generic phase diagram of the above model has 
been worked out in very recent publications^^ and is 
not repeated here. Our work concentrates on the metallic 
SDW and ACL states with small Fermi pockets, i.e. for 
magnetic fields H bigger than the critical field strength 
H c2 above which superconductivity is destroyed. 

In the above model we again introduce disorder and 
magnetic field on the level of the single particle propaga- 
tor. We realize that there are two cyclotron frequencies 
uj cg and tOcf, associated with the two different sorts of 
electrons. Additionally, there are two scattering times 
T g and Tf. The following discussion is again valid in the 
limit iv cg Tg,Lu c fTf <C 1. The specific heat of the system 
is this time composed of four different contributions 



C V 



N z C Z y 



Nc 



(51) 



The first term is due to the bosonic spinons and was cal- 
culated in Ref^. The two following terms were calcu- 
lated together with the magneto-oscillations in Eq. (|30| . 
Turning to the gauge field propagator we realize that the 
gauge field propagator is still given by Eq. 1(5)), however 
Eq. ([9]) modifies to 



IT 



mi 



v f 



Dfl 2 



D f q 2 



n yy 

which implies we can write 

b- 1 = D 



luoy + q Xg+ luay +q Xf 



ioja x g v + iuaf 



(52) 



(53) 



9 ' "f ■ 

A further modification comes into the picture due to the 
presence of the bosonic spinons. In the SDW state they 
condense, i.e. (z) 2 ~ Aaf, implying that the Higgs effect 
contributes a mass term causing U vv = iioa v g v + q 2 Xg + 
iuia v f v + q 2 Xf + A A f- 

For the moment, we will neglect this term and defer 
the discussion of the SDW state to the end of the sec- 
tion. As we discussed earlier, the oscillatory contribu- 
tions to the thermodynamic and transport quantities en- 
tering Eq. (152|) are of the form e ^<=s T s and e W . In 

the following we will consider e -C e ^"^a . With 

this we can express the oscillating term to leading order 
in — according to 



D~ 



2-KU, _ 

-2 cos e 



(i 



24 -^cos^e-*C(l> cs ) 





o x g q 2 



o 




(54) 



As in our preceeding discussion in Sec. IIVB 21 the first 
term drops out and eventually we can write the whole 
expression as 



Imtr D -D- l c 



-24- 



Cq b co + Cq 2 0J 3 



Aq s 



/i 2-ITfJ, - 

— =- cos e 

Aq^UJ 2 + BUJ 4 UJ 2 g LU r . 



J >:g 



(55) 



Again, we are only interested in isolating the logarithmic 
behavior, for which we only need to know the constants 
A, B, and C. We find following the earlier calculation 
that the logarithmic part of the specific heat is given by 



C n 2 2irfi _ 
—r — 5- cos e 

A ujig UJ r 



(56) 



with 



A = D 2 D 2 



vx 



B = {a vv Y , 
C = l)~l)ii' 2 " ! 



Xg 



(57) 



where 



V 

X 



Vg + Vf , 

Xg + Xf ■ 



(58) 



In the SDW phase, the presence of a finite Higgs term 
Aaf 7^ in the gauge field propagator introduces a new 
energy scale into the problem and the character of the 
problem changes its character to that of a Fermi liquid 
and the magneto oscillations will essentially be given by 
LK theory. Without going into the details we can show 
this along the lines of the derivation of Eq. (f56|) . The 
logarithmically in temperature diverging prefactor is now 
cut off by the Higgs mass A^- A compact formulation 
of Eq. (|56[) treating both regimes is given by 



C 1? 2irn _ 
— —5- cos e 

A uj£ g UJ r 



"eg 



x In 




A 2 



max 



T^T 2 /(T,., 



eg) 1 



(59) 

which in the SDW phase at very low temperatures (T <C 
A £ F ) corresponds to a simple renormalization of Fermi 
liquid theory. 
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VI. SUMMARY AND DISCUSSION 

Within this paper we analyzed magneto-oscillations 
in the specific heat of two-dimensional non-relativistic 
fermions coupled to a U(l) gauge field. Our results apply 
to a variety of different problems, including the descrip- 
tion of the 1/ — 1/2 fractional quantum Hall stat o 19 ' 26 i 27 
and different effective low energy gauge theory descrip- 
tions of the cuprate superconductor a 10 i 28 ' 29 i 30 . Our main 
result obtained in this paper is the calculation of the cor- 
rection to the standard Fermi liquid result of Eq. (|3T))) . 
which is shown in Eq. (|44f . We point out a couple of 
generic situations, in which our result holds and further- 
more discuss its meaning in the context of the gauge the- 
oretic description of the underdoped cuprates introduced 
by Galitski and Sachdev^. 
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APPENDIX A: THE GRAND POTENTIAL 

The grand potential of the disordered Fermi gas in a 
magnetic field is readily calculated using the fermionic 
propagator 



r^EM-c^K.^)) (ai) 



where 



G(ui n , E m ) 



(A2) 



with E m — u) c (m + 1/2). We use the following Poisson 
summation identity 7 - 

?'(-H))-f>> 



1=1 



with 
fix) 



, . 2nlx 
t dxj (x) sin 

ZTTl In UJ r 



-1 



(A3) 



lLU n + (J, - X + l^f- 



1 

2ttt 



du 



iu n - oj (oj + (i - x) 2 + -rK 



(A4) 



We can now perform the Matsubara sum yielding 

/CO 
dumj(u)G"(w) . (A5) 
-°° 

We furthermore use ttKt = vq.uj c , which implies 



Q = -/3 _1 ^!/ j dE \n(-iui n + fi - E 

n 

CO 

+ 2^ w c ^2 ■ 



. sgnw„ 
2r ' 



' dE sin x 



1=1 



2tt 2 / 



x / dujn f {uj)G"{uo,E) . 



(A6) 

We will now concentrate on calculating the second term. 
We integrate the energy integration by parts to obtain 



f 

Jo 



dE sin / dton f (Lu)G"(uj,E) 

u c J 

dunt(w)£l cos o) - G"(u, oo)) 

2nl UJ C 



LOr dE 



(A7) 



The first term on the right-hand only produces one term. 
This can easily be checked, since for E = the expres- 
sion is finite, whereas for E — > oo one can see, that the 
term vanishes, since the Fermi function restricts u> to be 
smaller than zero. Consequently, we obtain 



dE sin^^ / duon f (uj)G"{uj,E) 

,1 UJ C 



J dum/(u)|^G>,0) 



, s LUc f°° 2nlE d .,, 

dujn f (uj)-^- dEcos —G"(u,E). 

zttI J Q oj c dE 

(A8) 

We can split the remaining task into two parts. For the 
first integral one obtains in the limit /i 3> 5, T 



dw n / ( w )^G>,0) = -^ 



(A9) 



We now calculate the second integral 

OJ, 



dumf(uj) 



2ttZ 



o 



/•OO 

dnf(u>) lu 



duj- 



dw 2~kI 



2ttIE d .,, 

dEcos — G"(u,E) 

oj c dE 

cos 2^(G"( W , -oo) - G"(w,oo)) 
dEcos^^G"(uj,E) . (A10) 



The first term again is easily analyzed. For u> — > oo 
the Fermi distribution annihilates the expression. For 
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u> — > — oo we realize that the denominator of the Green's 
function diverges, since E > 0, which implies that the 
denominator overall goes like — : — h-. — j— for uj — > — oo, 

thus going quadratically to zero. This implies 



= / dio 



, ,u c f 00 , 2nlE d „„, 

dum f (u)-^ / dE cos —G"{u,E) 

2-kI J Q 

diif(uj) oj c '' v 



dio 2-nl 



dE cos ■ 



lu c dE 
2nlE 



G"(u),E) 



(All) 



One can easily obtain the -j- — > limit of this expression. 
This yields 

24 277^ oo c 

(eB) 2 (eB) 2 2iru „ N 
= fi + i_^-^-^cos— ^e"— . A16 
487rm Air^m uo c 

APPENDIX B: THE FREE ENERGY OF THE 
OSCILLATORY PART OF THE PHOTON 
SYSTEM 



This expression can be treated by noting the derivative of 
the Fermi energy pins the w-integral to zero. This allows 
to treat the energy integration according to 



In the following we sketch the isolation of the temper- 
ature dependent part of the oscillatory part of the free 
energy. We start with an expression of the type 



I) 

a dEcos 2nliE + ^G"(u,E-n) 



dEco S 27Tl(E + fJ ' ) G"(u J ,E-ri 



2nlfi 2txIlo 
—it | cos cos 

L0 C L0 C 



, 2ixl[i . 2ttIlo 
sm sin I e 



(A12) 



fo SC = a I 1 d< l duon b (uo) 
Jo J -A' 



uoq 



Aq 4 + to 2 



(Bl) 



where A w 2k f and A' w Ef. We first perform the 
integration with respect to q 



Zoic = J duon b {uj)uj\n 



( AA 4 + u 

and an integration by parts which leads to 



(B2) 



From there we can go on to solve the remaining integral. 
Since sin 1--- is an odd function of oo and d "f ^ even, 
the integral over sin drops out leaving us with 



dntiuj) lu c 2-kIii 2-kIuo 
dio — — cos cos e 



dio 21 oo c 
21 C ° S oo c 6 " T sin h2^Z 



(A13) 



We finally obtain (except for the diamagnetic contribu- 
tion this can be compared to Ref. 31 ) 

n n VoU) c (-!)' (-, 2lrl V X i —is 
il = il — TT^r y. — - — I 1 — cos , ; : e 



uo c sinh A; 



Oo 



2tt 2 I 2 
i=i 

v ujI vquj 2 c y> (-1) ' ^ 2nlfi Xi ^jj, 
24 2tt 2 ^ I 2 



■ cos • 



oo c sinh Xi 



(A14) 



where Xi = 2ff lT . In the limit oj c t < 1 we can approxi- 
mate the grand potential as 

n = n + — cos — — e <*ct . (A15) 

24 2tt 2 uj c sinhAi V ' 



A- WW* 



aA 4 
a 

8l 



In (AA 4 + (A' 

f A ' , /, 2, f, ^A 4 
duon b [oo)uj In 1 H ^~ 

-A' V w 



A 4 " A ' 



-A' 



dum^) In {lo 2 + ,4A 4 ) . (B3) 



The first two terms can be discarded, since they have 
no temperature dependence. The last term will also not 
contribute a temperature dependent part to leading order 
in £r and X . Proceeding with the remaining parts of the 
integral we obtain 

Josc 8A T J_ oc 1-coshw 

+ ^-T 2 r du;^^ 
8A /_„ 1 — coshw 



_£L ^ T 2 ln _ o.05880^-T 2 . (B4) 
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